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Abstract 



It is shown that at low densities, quantum dots with few electrons may 
be mapped onto effective charge-spin models for the low-energy eigenstates. 
This is justified by defining a lattice model based on a many-electron pocket- 
state basis in which electrons are localised near their classical ground-state 
positions. The equivalence to a single-band Hubbard model is then established 
leading to a charge-spin (t — J — V) model which for most geometries reduces 
to a spin (Heisenberg) model. The method is refined to include processes 
which involve cyclic rotations of a "ring" of neighboring electrons. This is 
achieved by introducing intermediate lattice points and the importance of 
ring processes relative to pair-exchange processes is investigated using high- 
order degenerate perturbation theory and the WKB approximation. The 
energy spectra are computed from the effective models for specific cases and 
compared with exact results and other approximation methods. 
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1 



Typeset using REVTeX 



I. INTRODUCTION 



Technological advances in microfabrication with corresponding reduction in feature sizes 
has led to renewed interest in the transport properties of semiconductor submicron structures 
both from the fundamental physics point of view and for possible future applications. In 
small metallic islands the Coulomb blockade gives rise to the single-electron effects, which 
may be modelled classically by a small intra-dot capacitance Q. Islands fabricated on 
the basis of semiconductors, called quantum dots, show in addition to the charging effects 
discrete energy levels related to size quantisation fH-|J. Individual quantum dots may be 
fabricated from hetero-structures, which limits the electron motion to two dimensions, by 
imposing lateral confinement with a metallic gate electrode deposited using fine-line (e-beam) 
lithography M. The discrete energy spectra may be measured in transport experiments 
at finite voltages or frequencies 0. The influence of Coulomb correlations on the 
excitation energies is difficult to see optically by far infrared absorption || due to the 
generalized version of Kohn's theorem ||. The observation of the quadrupole transitions by 
using grating couplers has been suggested |10| as a possible optical means. 

Quantum dots have been referred to as 'artificial atoms' [JTXJ] since the number of con- 
ducting electrons can be smaller than ten or twenty. At the low electron densities which 
are experimentally accessible, the calculation of excitation energies becomes a challenging 
theoretical problem. For real atoms, an independent-electron picture correctly describes the 
main physics and the Hartree Fock approximation yields reasonably accurate eigenstates 
which may be refined in a controlled way by perturbation theory. This is not generally the 
case for semiconducting quantum dots for which even an optimal Hartree Fock approxima- 
tion is significantly in error [12| and can even give qualitatively incorrect results, such as the 
wrong spin- mult iplet structure. The reason for this qualitatively different behaviour is that 
the electrons in a semiconducting quantum dot are highly correlated, due to the low effective 
density and the restriction of the electron motion to only two dimensions. Many body effects 
must be taken into account. This has been done by numerically exact diagonalizations for 
systems with a very few electrons (N < 4) ||13|| . In order to obtain reasonably accurate 
spectra for systems with more electrons, approximations must be made. For electrons in 
high magnetic fields, the low-energy states of maximal spin have been determined to good 
accuracy by making a harmonic expansion about classical minimum energy configurations 
for circular dots with parabolic confining potentials [[14]]. The resulting effective Hamiltonian 
is then diagonalised explicitly and electron antisymmetry is subsequently imposed to deter- 
mine the allowed eigenstates. High accuracy is achieved by including a sufficient number of 
Landau levels in the basis set. A complementary technique which has achieved some measure 
of success at low electron densities, and is not restricted to high magnetic fields, is based on 
many-electron "pocket" states JL5|. Low-energy spectra of systems with up to six electrons 
have been determined accurately |TB[ and the agreement with the exact numerical solutions 
(N < 4) is good. Unfortunately the computational effort is likely to become prohibitive for 
systems with more than ten electrons and we are again faced with the problem of devising 
a reliable approximation method for such cases. 

One possible route is suggested by the pocket-state analysis itself which exploits the per- 
mutation symmetry of the wavef unctions. We notice that in many cases the eigensolutions 
may be represented by an effective spin model of a simple form (Heisenberg model). This 
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apparent equivalence is reminiscent of magnetic insulators for which the correlated electron 
problem is known to reduce to a spin-Hamiltonian, a mapping which may be justified by 
transformation theory W% or degenerate perturbation theory [JTS | as well as, in some cases, 
by the theory of permutation groups ||19|| . The magnetic insulator problem starts with a 
lattice model in which the electrons are localised on atomic-like orbitals and for which elec- 
tron correlations are essential, the generic model being the so-called Hubbard model [20|. 
The underlying crystalline lattice establishes sites on which suitable one-electron states are 
centred. The question arises as to what circumstances, if any, might a Hubbard-type model 
be applicable to the problem of interacting electrons in a quantum dot, for which the Fermi 
wave lengths are much larger than the interatomic spacings and the underlying crystalline 
lattice loses its significance. At extremely low electron densities the long range interaction 
energy, according to Wigner JHJ], creates a crystalline like ground state which might define 
new electronic lattice sites. This crystallization is expected to take place only at electron 
densities which are a factor of 20 smaller than the densities used in experiments. Never- 
theless, it has been shown that the assumption of localised electrons allows the calculation 
of discrete low energy excitations as quantum corrections to the Wigner crystal energies 
fT6fl . The validity of a Hubbard model description is not only an interesting question of 
principle but has the practical potential of enabling systems of many more electrons to be 
dealt with than has hitherto been the case, i.e. tens of electrons or more. This would en- 
able the full machinery of techniques for solving the Hubbard and Heisenberg models to 
be applied, for which there has been much progress in recent years following the discovery 
of high-temperature superconductivity and the associated theoretical activity on correlated 



electron systems [22 



In this paper we will show that such a mapping may indeed be justified for the a priori 
continuous problem of interacting electrons in a quantum dot at low densities. This is 
done through the use of pocket-states which are briefly reviewed in the next section after 
introducing the basic interacting electron model. In section |TJ we show how the pocket- 
states may be used to define a one-electron orthonormal basis with orbitals localised on a 
"lattice", related to the electron configuration(s) in the classical ground state (determined by 
electrostatics). This is then used to construct a tight-binding Hamiltonian which reduces to 
the Hubbard model at low electron densities. Further mapping to an effective (Heisenberg) 
spin Hamiltonian or charge-spin (t — J — V) model is then performed. 

In section [IV] we point out an essential difference between quantum dots and the cor- 
responding correlated electron problem for interacting electrons on a true lattice of atoms, 
namely the increased importance of the so-called "ring" terms in the former case, which 
involve cyclic permutations of more than two electrons. The simple Hubbard model on a 
lattice underestimates the magnitude of these ring terms for the dot and it is shown how 
the method may be refined by introducing intermediate lattice points which are unoccupied 
in the ground manifold but give rise to the required ring processes through virtual excita- 
tions. In section [V| we solve the effective Hamiltonians for specific examples with N < 6 
and compare with exact numerical and pocket-state results where appropriate. Finally, in 
section [Vl|, we give a summary and discuss the outlook for this approach in dealing with 
systems of more electrons and of obtaining further corrections where necessary. 
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II. THE MODEL AND POCKET-STATE BASIS 



We consider the iV-electron quantum dot described by the Hamiltonian 

H = E (|| + v ( x ^ + W ( x i ■ ■ ■ x n) (2-1) 



W{x l ...x N ) = -Y J w{\x i -x 3 \) , (2.2) 

where Xi and p { are position and momentum of the % -th electron in d-dimensions (d = 2 
for most quantum dots) with (effective) mass m and spin s = 1/2. Neither the one-particle 
confinement potential v(x) nor the interaction w(x) depend explicitly on spin. 

The cases of square well single particle potentials in one dimension (ID) and in two 
dimensions (for square geometry) have been discussed in To be specific, and for com- 
parison later, we shall also mainly confine ourselves to the 2D square potential well of square 
geometry, though the extension to other geometries is straightforward and, for the present 
consideration, the detailed form of v(x) is not qualitatively important. At large mean 
inter-particle distances, r s , it becomes energetically favourable for the electron system to 
localize its charge density distribution in regions close to the classical ground-state electron 
configuration(s) [plp3| , which may be determined by minimizing the electrostatic energy. 

This fact is the motivation for the many-electron pocket-state basis in which the Hilbert 
space is restricted to 1 < p < v ■ N\ basis states \p) with spatial representations being 
defined in configuration space of dimensionality d N . The number of classical minimum 
energy configurations v may be greater than one in certain symmetric geometries. (For 
example, for N = 3 on a square we have v = 4 since, for the classical ground-state energy, 
any one of the four corners may be unoccupied.) The low-energy levels form a multiplet 
which contains v ■ 2 N states in total (including Zeeman levels). The excited 'vibrational' 
states scale like a power law ~ r^ 7 with r s ( 7 being close to 3/2 ) and can be neglected 
if the density is not too high. This is due to the scaling ~ exp(— \Jr s /r c ) of the low 
energy excitations we are interested in here, where r c is a length scale which characterizes 
the transition from the almost non-interacting situation (r s <C r c ) into the qualitatively 
different regime of strong correlations (r s 3> r c ) ]n| . 

In the pocket-state method, the first step is to calculate energies by assuming equivalent 
but distinguishable particles and ignoring their spin. Spin and statistics are subsequently 
reestablished by means of group theory. Tunneling integrals (p\H\p') between the localized 
many particle (pocket-) basis states determine the low-energy excitations. The tunneling 
processes correspond to (correlated) transitions between different particle arrangements. In 
the low-density regime one tunneling integral, J/2 , is exponentially larger than all others. 
In many cases the dominant J corresponds to the exchange of only two particles, which 
may be adjacent in real space. Other processes could be the simultaneous (ring-) exchange 
of three or more particles, corresponding to a cyclic permutation. 

An important feature of the spectra given by the pocket-state approximation (and by the 
effective charge-spin models to be derived and discussed in the following sections) is that 
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the fine-structure of the energy spectra depends only on J, the magnitude of which may 
be estimated semiclassically within the multidimensional WKB approximation. The ratios 
between the energy differences are insensitive to the detailed form of the inter-electron po- 
tential and to r s . An examination of the spectra shows that these situations may be mapped 
onto an effective spin model with again one parameter for each process (pair exchange, ring) 
considered. This equivalence will be justified in the next section. Here we merely make the 
observation that for the cases considered by the pocket-state method, the energy spectrum 
is the same as that which would be given by a spin-Hamiltonian. An exception to this are 
situations for which more than one classical minimum energy configuration for the electrons 
exist ( v > 1 ), such as are 2 or 3 electrons on a square. We will show in the next section 
that these latter situations may be mapped onto a so-called t — J — V model in which the 
dominant tunneling process is now related to t, the amplitude corresponding to the jump of 
a single particle into an empty 'site'. J is again an exchange process involving at least two 
particles and V represents the Coulomb repulsion between electrons on neighboring 'sites'. 

In general, the pocket-states are correlated electron states which need not be assumed as a 
direct product of one-electron states. However, we may mimic the eigenstates at sufficiently 
large r s with a one-electron product basis within the Heitler London (HL) approximation. 
The latter are localised in regions close to the positions the electrons would have in the 
classical ground state. In the next section we show how this one-electron basis may be used 
to define an orthonormal basis from which we construct an effective Hubbard model. This 
effective Hamiltonian yields low-energy spectra which are at least as accurate as what would 
be given by the corresponding HL approximation using a product basis of non-orthogonal, 
one-electron wavef unctions. We then show in section [TV] how the method may be extended 
to cases where ring processes are important. 

III. MAPPING TO HUBBARD AND CHARGE-SPIN MODELS 

Let us assume that the pocket basis states \p) are approximated "optimally" by non- 
orthogonal one-electron wavefunctions 0, centred at position i. These states have the form, 



where {p 1; . . . ,Pn} is the permutation p of the sequence {1, . . . , N} . In this section we 
orthogonalise these one-electron wavefunctions and construct antisymmetrised A-electron 
states from which we define the effective Hubbard model, which is subsequently transformed 
into a spin-Hamiltonian. To illustrate the main features of the method with minimal mathe- 
matical complexity, we first consider in detail the simplest non-trivial case of two electrons in 
a ID square well. The analysis is then generalised to A-electrons in one or two dimensions. 

As is well known, for two electrons the antisymmetrised states may be written as the 
product of orbital and spin parts, as with the HL states of the hydrogen molecule. The 
symmetric (singlet) and antisymmetric (triplet) orbital states are thus, 




(3.1) 



*5 




(3.2) 
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and 



1 (1)0 2 ( 2)-0 2 (1)0 1 (2) 
2(1 -s 2 ) 



(3.3) 



where s = (0i|0 2 ) is the overlap and the arguments (i = 1,2) abbreviate the coordinate 
Xi of the 2-th particle. These HL states yield approximate singlet and triplet ground-state 
energies: -Esingiet = (tys\H\^s) and triplet = {&a\H \^a), where H is the Hamiltonian for 



two electrons given by equation ( 2.1 



We now transform to orthonormal one-electron states ip where 



+ 



1 



2x/TTs 2vT= 



+ 



2 V / TT^ 2^/l^~S 



(3.4) 



'02 



\2^/T+s 2y/T=aJ T \2y/T+S 2^Y^S ) 
Inverting this transformation and substituting into equations (|3.2|) and 



we get, 



(3.5) 



(2) 



(3.6) 



and 



(1 ) ^(1)^(2) + ^(1M(2) 

v/2 

(2) ^(1)^(2) + ^(1)^(2) 

V2 



^(1)^(2)- ^(1)^1(2) 



(3.7) 
(3.8) 

(3.9) 



We note that ^ (triplet) has the same form as for the non-orthogonal states given in 
equation (|3.3| ), i.e. it corresponds to one electron in each orbital. On the other hand ty$ 

(singlet) consists of two components, one with an electron on each orbital ) and the 

(2) 

other with both electrons in the same orbital (W s ). (The latter state is not allowed for 
the triplet by the Pauli principle.) This "double occupation" is thus a direct consequence 
of orthogonalisation. 

Using these two-electron states composed of orthogonal one-electron wavefunctions will, 
of course, yield exactly the same approximate singlet and triplet energy expectation values. 
However, we can obtain a more accurate estimate for the ground-state energy of the singlet 
by diagonalising the Hamiltonian in the Hilbert space defined by the two-electron base 
states \E , 5' ) and (Eqn. |3]7]). This follows from the variational principle which ensures 
that diagonalisation of the Hamiltonian matrix will yield the optimum superposition of the 
base states, whereas the HL state ( |3.6| ) is not optimum in general. 
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This principle also applies to the general iV-electron stated in the introduc- 

tion, i.e. diagonalisation of the Hamiltonian in the restricted Hilbert space defined by base 
states (Slater determinants) constructed from orthogonalised one-electron wavefunctions 
will yield a more accurate low-energy spectrum than the simple products of non-orthogonal 
one-electron states ( |3.2| ) and (|3.3|). 

For the two-electron case the (2 x 2) Hamiltonian singlet matrix is, 



2s + V+j 
2t* 



2t 

2e + U + j 



where 



v = J \Mx)\ 2 \My)\M\x-y\)a 3 xd 3 y, 

U = J \Mx)\ 2 \My)\M\x-y\)d 3 x d 3 y = J \Mx)\ 2 \^ 2 (y)\ 2 w(\x-y\)d 3 xd 3 y, 



t= Ul(x) 



PL 

2m 



+ v(x) 



ip 2 {x) d 3 x, 



r = / 



\My)\ 2 w(\x - y\)d 3 y 



ip 2 {x) d 3 x, 



t = t + T, 



j = J ipl(x)ip;(y)w(\x - y\)ifo(x)tl) x (v) d 3 a; d 3 y, 

i = j ^l(xW 2 (y)w{\x - y|)V>i(a;)V>i(l/) d 3 x d 3 y, 

J=Re(£) . 

Note, j = £ = J if the tpi are real. 
The triplet energy is, 

K = 2s + V . 



3.10) 

3.11) 
3.12) 
3.13) 

3.14) 

3.15) 
3.16) 
3.17) 

3.18) 
3.19) 



(3.20) 



Within this restricted manifold of states, the singlet-triplet Hamiltonian matrix, equa- 
tions ( 3.10|) and ( |3.20|) , are equivalent to an effective Hamiltonian, 



H cS = e [m + n 2 ] + U [n n n n + n 2r n 2l ] + Vn x n 2 c L c laC 2a c m 
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+ 



t cLcfcr + T ( c Lc 2 a n ^ + cLc la n 2 a) + -f c\ a c\- a c 2 - a c 2 ^ + he 



(3.21) 



where c, 



1,2 satisfying (x\c\ a \v&c) = ipi a (x) 



is a Fermi creation operator at 'site' i 
n ia = c\ a c ia and = n iT + n^. 

This effective Hamiltonian is very similar to that considered by Hubbard for a pe- 
riodic array of one-electron atoms in the study of the metal-insulator (Mott) transition. 
Indeed, for this two-electron case, equation ( |3.21| ) is an effective Hamiltonian for the hydro- 
gen molecule within a restricted Hilbert space of (effective) ls-orbitals. We stress, however, 
that in the quantum dot case the localised orbitals are fundamentally different from those 
in the atomic case in that their very existence depends on the electron-electron repulsion 
which is responsible for localising the electrons near specific points in real space at low den- 
sity. For this reason we are not justified, a priori, in dropping all but the largest Coulomb 
term (U), as is usually done in the Hubbard model. In particular, the nearest-neighbor 
effective Coulomb interaction, V, becomes most important when r s becomes very large 
and equation ( |3.21|) simply reduces to the classical expression for the ground-state energy, 
independent of spin, as it should. This may be seen more clearly by transforming (|3.21 ) 
into an effective spin-Hamiltonian which is demonstrated most directly for this two-electron 
case by diagonalising the singlet-matrix, ( |3.1(J| ), explicitly to yield eigenenergies: 



E. 



± 



As + V + U + j + j ± J(U - V + j - j) 2 + 16\t\ 2 



(3.22) 



Since for the cases of interest (large separation between the electrons and real one-particle 
wave functions), U ^> V,j,3 and j = j , then, 



E 4 



2e + U + 3 + 



A\t\ 



U-V 



and 



2e + V +3 



A\t\ 



U-V 



(3.23) 



(3.24) 



Combining these equations with (|3.20| ) for the triplet we see that the energy spectrum 
consists of a singlet-triplet pair at low energies separated by a singlet at energy ~ U 
higher. Hence the low-energy singlet-triplet is equivalent to a spin system with effective 
spin-Hamiltonian, 



H, 



spin 



2e + V + J( Sl - s 2 -l/A) 



(3.25) 



where 



J 



A\t\ 



U-V 



u 



-j 



(3.26) 



The first term in equation Q3.26|) , favouring a singlet ground-state, is sometimes referred 
to as 'super exchange' and is usually larger than the second 'direct' exchange term, j. This 
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superexchange contribution has its origins in the assumed non-orthogonality of the initial 
Heitler-London basis (pi in (3.1). From the pocket states it is known that they form a 
non-orthogonal basis set for the low energy eigenstates |24j . 

We now consider the general case of N electrons in a quantum dot. Starting with a 
suitable product-basis of pocket-states, equation (|3.1|), we may, in general, define an or- 
thonormal basis by Lowdin's method [25). (For high-symmetry situations it may be more 
convenient and expedient to use some other method. For example, in the case of four elec- 
trons on a square we have 7r/4 rotational symmetry and we may generate orthonormal states 
by forming molecular orbitals, normalising and then transforming back to localised states. 
Similarly, for translationally invariant systems we may transform to running (Bloch) waves, 
normalise and then transform back to localised (Wannier) states.) Within the correspond- 
ing many-electron Hilbert space we may express the Hamiltonian in second quantised form. 
We are led directly to an effective 'single-band' model of the type considered by Hubbard, 
generalising equation (|3.21|), 



i j a 



+ 



13 o 



(3.27) 



where the indices i and j should be distinct. The parameters e, U, V,j, t, r and £ have the 
same meaning as in the two-electron case ( [3.11D , though it should be noted that they will, 
in general, depend on the position of the electron, or pair of electrons and on separation. 
In particular the parameters j, t, r and £ related to electron transfers decay rapidly as the 
distance between the 'sites' % and j increases. This allows us to retain only nearest neighbor 
terms in the double summations. As with the case of two electrons, we may eliminate the 
high-energy states of H c g to yield an effective spin-Hamiltonian. This may be done using 
degenerate perturbation theory |26| or a canonical transformation [27|, as with the case of 
a lattice. 

Specifically in the situations of two (or higher) dimensional quantum dots the cases of 
more than one minimum of the classical electrostatic energy may occur ( v > 1 ). Then the 
quantum mechanical charge density distribution shows more than N peaks at large r s , 
say N > N . The number of single electron states must be equal to N in these cases. 
N — N lattice sites will be unoccupied at an instant. Consider, for example, the case of a 
square confining potential in two dimensions. For N = 2, N = 4 and the electrons will be 
found mainly near opposite vertices along the diagonal of the square, due to the Coulomb 
repulsion. States with two electrons near adjacent vertices will be of order V higher in 
energy (nearest-neighbor Coulomb repulsion energy) and states with two electrons located 
near the same vertex will be of order U higher in energy. Similarly for N = 3 there will be 
low-lying states with electrons located near three of the N = 4 vertices with states having 
two electrons near the same vertex again being ~ U higher in energy. For iV = 4 all vertices 
are occupied in the ground manifold, with configurations for which there are two electrons 
near one vertex being at least ~ U higher in energy. 

In each of these cases the high-energy states, corresponding to one or more vertices 
being occupied by two electrons, may be eliminated from the Hilbert space by degenerate 
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perturbation theory, where they appear as intermediate states. More precisely, we first write 
H eS in the form H eS = H + Hi where H = [ti n i + Uin^riii}. H is diagonal in the basis 

i 

of product states with the highest-energy states having two electrons in the same state, 
being 



U higher than the low-energy states 
results in an effective H — J — V Hamiltonian 



Eliminating these states to second-order 



TliTlj 



P 



(3.28) 



In this equation t and V have the same meaning as before and 



J, 



<J 



2 \tt%j H~~ T{' 



1 



i 



Jij 



(3.29) 



P is a projection operator which eliminates all base states for which there are two electrons 
in the same localised state. (These states appear as intermediate states and give rise to 
the 'superexchange' in equation (|3.29|) .) We also note that the e, r and £ terms in equation 
( |3.27| ) have disappeared in equation (|3.28|) . The former has been dropped since it always 
gives rise to a constant and the latter two are precluded by P (as is the [/-term) since they 
would always give rise to double occupation at a 'site'. 

Equations ( |3.27|) and fl3.28j ) are the fundamental effective Hamiltonians for interacting 
electrons at low-density in a quantum dot and capture the main physics of the problem for 
most situations. (They do not, however, include processes in which groups of electrons can 
rotate simultaneously, which can be important in some circumstances. The generalisation 
to include these so-called ring terms is given in the next section.) Although these effective 
models were justified by starting with the Heitler-London approximation to pocket states 
using non-orthogonal one-electron states, we emphasize that they have a greater range of 
validity than the initial approximation suggests. The above methodology enabled us to de- 
fine a basis set of product states (Slater determinants) of some one-electron states and this 
led directly to effective Hamiltonians operating within this restricted (incomplete) basis set. 
However, it is well known , that higher- lying base states may be accounted for by pertur- 
bation theory. This technique has recently been applied very successfully to lattice models 
in which an effective single-band Hubbard model was derived from a multi-band model using 
quasi-degenerate perturbation theory pOfl . We shall not attempt such a reduction here for 
the quantum dot problem, but merely point out that its main effect would be to renormalise 
the effective parameters (e^, U^ty etc) and extend the validity of the resulting models. 

Let us apply (|3.28| ) to the case of electrons in a square well in two dimensions, discussed 
above. For N = 2 and N = 3 the t-term will cause electrons to hop from occupied to 
unoccupied sites whereas the J and V terms will be effective only when two adjacent sites 
are occupied. The presence of the V^-term ensures that for N = 2 the electrons will have 
lower energy when they are located on diagonally opposite vertices of the well, as mentioned 
earlier. For N = 4 the i-term may be omitted since all four vertices are occupied in the 
ground manifold and it would thus give rise to double occupancy, precluded by P. The 
V-term may also be dropped since it gives rise to a constant contribution. Hence, the sites 
i and j can be restricted to nearest neighbors so that ( 3.28|) reduces to a Heisenberg model, 
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Hj = J ij S i ' S j ( 3 - 30 ) 

<ij> 

where P has been omitted for brevity. 

Similar considerations apply for N > 4 and the effective Hamiltonian will either have 
the form (|3.30|) or ( |3.28|) , depending on whether v = 1 (e.g. N = 4, 5) or v > 1 (e.g. 
N = 2,3,6). For the latter further reduction will be possible when some of the 

states of the t — J — V model are clearly higher in energy due to the ^-interaction. 

Consider, for example, the case two electrons in a square-shaped dot for which there are 
six ways of distributing the two (indistinguishable) electrons in the four corners, resulting 
in a Hilbert space of dimension 4 x 6 = 24 when spin is taken into account. Four of 
these configurations (16 with spin) are of order V higher in energy than the remaining 
two, as shown in figure [I]. These high-energy base states may be eliminated from the 
problem by a further application of degenerate perturbation theory, resulting in an effective 
Hamiltonian which only operates in the Hilbert space of the 8 lowest-energy base states. 
By a straightforward calculation this takes the form (apart from an unimportant overall 
constant): 

H cS = A(nin 3 - n 2 n i ){R K/2 - R-^/2) + J[(si ■ s 3 - l/4)mn 3 + (s 2 • s 4 - l/4)n 2 n 4 ] (3.31) 

where A = 2t 2 /V and Rg is a rotation operator, i.e. 

R7T/2W, *, <?', *) = I*, o-, *, </) , (3.32) 

etc. Note that the prefactor (riin 3 — n-in^) takes the value ±1 depending on whether either 
'sites' 1 and 3, or 2 and 4 are occupied. This effective Hamiltonian is easily diagonalised. 
The spectrum consists of two singlets and two triplets. We see immediately that (|3.31|) 



gives zero when operating on a spin-polarised state and hence the triplets are at zero energy 
with the singlets at energies — J ± 2 A, in agreement with what is obtained from the pocket 
state approximation for this problem when the exchange of the two electrons is included 



there []16 |. There are thus two contributions to the 'binding energy' of the singlet ground 
state, a resonance energy (— 2A) and a superexchange energy (— J). It is easy to see that 
the former will be dominant since it occurs in second-order and does not involve double- 
occupation of a 'site'. On the other hand, there are two contributions to J, the second-order 
term ~ At\^/U and fourth-order terms ~ t\ 2 /V 2 U . We may extend these considerations 
to a regular polygon with N vertices and N/2 electrons. The effective Hamiltonian will 
again consist of rotation operators (of ±27r/iV) and exchange terms. The exchange terms 
are independent of N, whereas the rotation terms becomes relatively less important since 
they first occur in Nth-order. In the next section we shall encounter circumstances where 
simultaneous rotations of more than two electrons may be important. 

We conclude this section by emphasizing the computational advantage in this reduction 
to a spin or charge-spin Hamiltonian. Thus, for example, the Heisenberg model lies in a 
Hilbert space of dimension 2 N for iV electrons, which may be further reduced by exploit- 
ing other symmetries, e.g. conservation of total S z which reduces the largest subspace to 
-^V((t)') 2 (Sz = for even). This should be compared with 



max 

s 



(2S + 1)N\ 



(N/2 + S + l)\(N/2- S)\ 
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for the largest subspace of the Hamiltonian in the original pocket-state basis. 

The mapping also gives interesting insight into the nature of the energy spectrum. For 
those cases which reduce to a Heisenberg model, the spectrum consists of spin-mult iplets 
with energy-separation of order J ('bandwidth' ~ 2NJ). On the other hand, for cases which 
reduce to a charge-spin model, there are multiplets separated by energies of order t (>> J in 
the same geometry and electron densities), each of which have a fine structure with energy 
splittings of order J. 

The forgoing analysis for electrons may also be extended to a fictitious system of spin-1/2 
charged bosons. For repulsive Coulomb interactions they also show Wigner crystallization 
and the pocket-state method can be applied. This problem was considered by Hausler [16|] in 



order to proof the state of highest spin within the multiplet structure of electrons (fermions) 
in a quantum dot being the spin polarized state for v — 1 by exploiting an isomorphism 
between fermions and bosons using permutation group theory. We reconsider this problem 
here to enable another direct comparison with the pocket-state description. 

The case of spin-1/2 bosons is similar to the fermion case in that at low densities the 
largest 'Coulomb' matrix elements correspond to two bosons on the same 'site'. Neglecting 
all other matrix elements thus gives the 'boson' Hubbard model: 

H= £ + H.c.) + V - £ btblKK (3-33) 

<ij>a i,aa' 

Using the usual commutation rules for the Bose-operators yields 

b ] \n) = y/n + l\n + 1) 
b\n) = y/n\n - 1) . (3.34) 

It follows that: 

h \a, a') = U\a, a') 

where 

aa' 

This is the same result as the case of fermions except that with bosons the two particles 
may have the same spin. 

As with the fermion case, we can now reduce ( |3.33|) to a Heisenberg or t — J — V model 
(depending on whether or not the number of bosons equals the number of lattice sites). Let 
us, for simplicity, do this for just two sites. The extension to the general case may be done 
by perturbation theory in the same way as with fermions. 

In the 'atomic limit' (t = 0) the eigenstates of (|3.33 ) with two bosons consist of either one 



boson on each site with either spin (energy 0) or one site unoccupied and the other doubly 
occupied, again with either spin (energy = U). Using these base states it is straightforward 
to show that the singlet state: 

^ = — w~ 
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is unaffected by the hopping, i.e. 



H\il> s ) = 

On the other hand the triplet states mix. For example, using ( |3.34 ), we get: 

H\tv = t[\-,n) + \nr)} , 



tf|TT,-> = 0lTT,-) + >/2*|T,T) 

etc. This leads directly to the triplet Hamiltonian matrix: 



2t 
2t U 



with lowest eigenvalue: 



E 



U-VU 2 + 16t 2 



This is exactly the same result as for the fermions except that the positions of the (low- 
est) singlet and triplet levels are reversed. It follows that the lowest singlet and triplet is 
equivalent to a Heisenberg model with ferromagnetic exchange of the same magnitude as the 
case of fermions. As mentioned above, this may be generalised to any number of sites by 
degenerate perturbation theory resulting in a Heisenberg model in which J — > — J when 
compared with the fermion case, leading to an 'inverted' spectrum in agreement with the 
permutation group analysis described in the second reference of |Tj| . 

In a similar fashion, cases for which the number of bosons is less than the number of 
sites (such as 2 or 3 on a square) may be reduced to a t — J — V model. However, whilst 
the sign of J is reversed relative to the case of fermions, the sign of t remains the same and 
the resulting spectrum is no longer a simple inversion. 



IV. RING TERMS 

In the section [Hj it was shown that the t — J — V model may be reduced to an effective 
Hamiltonian which only involves the lowest-energy base states. This effective Hamiltonian 
then contains operators which rotate some or all of the electrons simultaneously (cf. (|3.31|) ). 
The possibility of such a collective motion of a subset or all of the electrons is much more 
general and will in fact occur for any geometry. These processes are negligible within the 
framework of the usual lattice models but this is not always the case in quantum dots or 
other semicondutor based nano-structures where such 'ring' terms can become dominant. 
As we shall see they may be included in a lattice model by introducing intermediate lattice 
points which are unoccupied in the ground manifold. 

Consider, for example, the case of an equilateral triangular dot containing three electrons, 
as shown in figure ^. Using the Hilbert space defined by the Heitler-London approximation 
to the pocket-states, leads directly to the Heisenberg model for this system, with a spin- 1/2 
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ground state, as shown in section [II. Ring terms corresponding to cyclic rotations of all 
three electrons are not included in this model. However, it is clear from the pocket-state 
description that such processes exist. The relevant tunneling integrals between the 3! = 6 
pocket-states in this case can be represented by permutations of the types 



(123) — ► (213) 



J/2 Pair Exchange 



(4.1) 



(123) — ► (312) 



K/2 Ring Exchange 



(4.2) 



Their magnitudes may be estimated semiclassically within the WKB approximation. The 
most important contributions to J and K vary exponentially with the particle distance 
r s (which equals the length of one side of the triangle in the case considered), i.e. 



J/2 



~ e 



Sn 



(123)-f(213) 



where 



K/2 



~ e 



-<S(123)-» 



(312) 



s 



dq J2mnW(x(q)) 



is the action associated with the transport of n particles (n = 2 or n = 3 for J or K , 
respectively) of mass m from the initial classical ground state (123) into the final state. 
The transition takes place along the trajectory x(q) which in principle obeys an equation 
of motion which minimises the classical action. Here, it is parametrized by < q < 1 in 
configuration space (of dimensionality dN = 6 ). For an explicit estimate let us specify the 

e 2 Ik 

interaction between the particles and assume w(\xi — Xj\) = -, r . Then W(x) is 

| *^ ' i *^ J I 



l «U^' *jC rj 

the Coulomb energy at x = {xi,X2,x$} where the equilibrium value 
subtracted so that W(x(0)) = W(x(l)) = . The results 



, has been 



S 



(123)-«-(213) 



4y m — r s 
V K 



1/2 



dq 



1 



1 



y/1 - q + q 2 



2'.) 



q + 4g 2 



7^2 



1/2 



(4.3) 



2.139Jr s me 2 /K 



S 



(123)-«-(312) 



bum — r s / dq 
V K Jo 



y/l-3q + 3q 2 



1/2 



(4.4) 
(4.5) 



= 2.852 v /r s me 2 //« (4.6) 

are obtained assuming straight lines x(q) for simplicity. This yields upper bounds to the 
true values of S which may be reasonable particularly for <S(i23)-+(3i2) • The estimate for 
5(123)^(213) is surely worse because the classical path x(q) is more difficult to estimate for 
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the pair exchange. The chosen paths for both cases are indicated in figure |2j. Thus the pair 
exchange process (corresponding to the superexchange within the lattice description below) 
is slightly dominant compared to the ring process in the equilateral triangle. This supports 
a ground state of low total spin S = 1/2 . However, the triangle can easily be distorted so 
that 5(123)^(312) is much reduced while <S(i23)->-(2i3) remains almost unaffected, resulting in 
a crossover to a spin polarised ground state |31~l , [32f (see also below). 

The existence of ring exchange processes highlights a fundamental difference between a 
true lattice model and a continuum model. In order to accommodate such processes we 
need a better lattice approximation to the continuum than that suggested by the Wigner 
lattice. This may be achieved by introducing intermediate lattice points and corresponding 
localised one-electron orbitals. Orthogonalising these one-electron orbitals then enables an 
orthonormal many-electron basis set to be constructed, as before. This is shown in figure |3| 
for the triangular dot in which an extra lattice point is inserted between each pair of vertices. 

The low-energy eigenstates may now be expanded in this basis set. As before, the prob- 
ability of occupation of the three vertices is highest, leading to a 2 3 = 8-fold degenerate 
ground-manifold when the hopping terms are switched off. Higher-lying base states, cor- 
responding to double-occupation of a lattice site and to occupation of intermediate lattice 
sites may then be eliminated by degenerate perturbation theory, resulting once more in an 
effective spin-Hamiltonian. Intermediate states which involve double-occupation of a lat- 
tice site give rise to a Heisenberg superexchange term in the effective Hamiltonian, with 
J = 8t 4 /V 2 U, where t, V and U have the same meanings as before but now refer to nearest- 
neighbor and on-site interactions for the new lattice. Longer-range interactions are ignored. 
(Note, however, that this is not necessary, it being straightforward to include longer-range 
interactions leading to a renormalisation of J which will leave the form of the (Heisenberg) 
pair exchange unchanged.) There are two processes involved in this superexchange making 
equal contributions to J. These are shown in figure |j. These processes are analogous to 
superexchange in atomic systems in which the intermediate sites would correspond to lig- 
and ions surrounding transition-metal ions. Recent examples of this kind of system are the 
copper-oxide planes of the high-temperature superconductors, for which there is a similar 
fourth-order expression for the superexchange |33| . It is also clear from this new lattice de- 
scription that other processes, which do not involve double occupation of a site (and therefore 
do not involve the exclusion principle), are potentially important. The lowest-order process 
of this type is also a pair exchange and occurs in fifth-order. An example is shown in fig- 
ure ^ (b), which makes a contribution to J of order t 5 /V A . Despite being higher-order than 
the other superexchange terms, this process is, in fact, the largest since it does not involve 
double-occupation of a site (i.e. there is no [/-factor in the energy denominator). 

With intermediate lattice points it is easy to see that there are ring processes giving 
rise to simultaneous rotation of all three electrons. One such process, in which the three 
electrons are rotated cyclically by 2tt/3, is shown in figure |5]. The process resembles the 
classical path which would contribute to the WKB tunneling rate, cf. figure |[ The former 
is a sixth-order process giving rise to a cyclic rotation operator in the effective Hamiltonian 
with amplitude t 6 /12V 5 . Summing over all such processes results directly in an effective 
'ring' operator, —K(R 27T /3 + R-2n/z)i where K is a constant energy and #2^/3 (-R-271-/3) * s an 
operator which performs a cyclic (anticyclic) rotation of all three electrons (cf. ( |3.32|) ). We 
note that these rotation operators may be written as spin-operators using the identities, 
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£W/ 3 = P12P23 and P_ 27r/3 = P 23 Pi2 (4.7) 

where Py = 2si • s 2 + 1/2 is the Dirac exchange operator. Hence, the effective spin- 
Hamiltonian for 3 electrons in an equilateral triangular potential well becomes, 

H cS = J - [(P 12 - 1) + (P 23 - 1) + (P31 - 1)] - ^(Pi 2 P 23 + P23P12) • (4.8) 

We see immediately from this form that the spin polarised state has eigenenergy —2K. 
It is straightforward to complete the diagonalisation of H e g which results in a further pair 
of degenerate doublets at energy K — 3 J/2. Thus we see that the ring and exchange terms 
oppose each other with the former favouring a high-spin ground state and the latter a low- 
spin ground state. The crossover occurs at J = 2K, which is not quite reached for the 
equilateral triangle, as discussed above. 

This will not always be so for other geometries as follows immediately from the fact that 
we can change the shape of the dot in such a way that the ring processes are increased in 
amplitude relative to the exchange processes. An 'extreme' example of this is the circular 
dot for which there is no barrier to the simultaneous rotation of the three electrons and 
the ground state should be spin polarised. The transition from high-spin to low-spin ground 
states in rings with increasing impurity barrier has been discussed recently |J2| in the context 
of persistent currents. 

With the inclusion of higher-order exchange processes involving also more than two 



particles, this completes the mapping of formerly continuous electron problems ( |2.1| ) onto 
lattice models fl3.28|) for their low energy and spin properties. This incorporates all per- 
mutational processes which according to the pocket-state description may be relevant. The 
above derivation of an effective-spin model for a triangular dot, including ring terms, is easily 
generalised to arbitrary shapes by inspecting all possible classes of ring and exchange terms 



on a suitably dense lattice into (]3.28|) . The relative amplitudes of the various processes can 
be estimated semiclassically. 



V. COMPARISON WITH EXACT RESULTS AND OTHER APPROXIMATIONS 



To corroborate the lattice description for the low energy and spin properties of quan- 
tum dots we compare the results from ( |3.28| ) with results obtained by numerically exact 
diagonalisations [[L5[] and by the pocket-state method . 



The simplest model for an "artificial atom" is a hard wall box in one dimension. Rigorous 
properties for the sequence of spin states 



E(S) < E{S') 



(5.1) 



for S < S' are known by a theorem of Lieb and Mattis [34]. E(S) denotes the lowest 
energy of N almost arbitrarily interacting electrons to given spin S = j ^2 j ^ • • • > N/2 . In 
particular, the ground-state spin must be minimal. This rule ( |5.1| ) is in obvious agreement 



with the description in terms of an antiferromagnetic Heisenberg (spin-1/2 ) chain ( [3.30D , 
valid for this case according to section [TIT] , where the exchange of adjacent electrons is the 
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relevant permutational process and there is only one classical minimum energy configuration, 
( v = 1 ). In table | the energy values of N = 2, 3, 4 electrons are given in units of J as 
obtained within the pocket-state approximation (PSA), assuming J 12 = J23 = J34 ■ These 
spectra agree exactly with the eigenvalues of the Heisenberg chain of iV = 2,3,4 spins 
apart from unimportant additive constants. Similarly, figure ^| shows the excitation energies 
for N = 3, . . . , 6 and includes also results obtained by numerical diagonalisation of the 
Hamiltonian matrices for N = 3,4 (dashed). For the cases N = 5 and iV = 6 we also 
checked, by numerical diagonalisation, that the Heisenberg chain shows exactly the same 
spectrum. 

As a further example we investigated N = 2, . . . , 5 electrons in a two-dimensional square 
with hard walls. For these cases the interaction must be long range w(|£c|) ~ \x\ ~ 7 with 
7 < 2 in order to make the PSA applicable. The spectra obtained within PSA are shown 
in figure [7]. For N = 4 the nearest neighbor pair exchange process has been compared 
with the competing process of all 4 electrons rotating cyclically by tt/4 flTEfl , in a similar 
way as it has been done for N = 3 electrons in a triangle in the previous section, with 
the result of the former being slightly dominant. In all cases the equilibrium positions of 
the electrons have been taken as sites for single electron states for the lattice description, 
as discussed in section [TTT| , and the importance of transition processes has been estimated 
semiclassically. For N = 5 one electron is, by electrostatics, located in the centre of the 
square being coupled antiferromagnetically, according to (|3.30|) , to the 4 electrons placed in 
the corners. The Heisenberg model is appropriate since there is only one classical equilibrium 
configuration ( v = 1 ). Neglecting next nearest neighbor exchange results in a S = 3/2 
ground state spin, using the PSA and permutation group theory. This result is immediately 
understandable from the Heisenberg spin-Hamiltonian, which clearly has a quartet ground 
state. Furthermore, we have checked that for all spectra, the corresponding Heisenberg 
( N — 4, 5 ) or t — J — V (N = 2,3) Hamiltonians exactly reproduce the correct ordering of 
the spin multiplets and the correct energy splittings provided that the parameters t, J and 
V are chosen to be consistent with the corresponding tunnelling amplitudes in the PSA. 

As a third example we look at the spectrum of three electrons in a one-dimensional ring. 
In section |TV] we explained why the ground state is high spin ( S = 3/2 ), which is also 
found in the pocket-state description. Two degenerate low spin ( S = 1/2 ) excited states 
are higher in energy by the rotational constant h 2 /2mNL 2 , where L is the circumference 
of the ring. They are additionally shifted by the magnitude of the pair exchange process J 
towards lower energy, in agreement with 



VI. SUMMARY AND OUTLOOK 

Based on the pocket-state description of the low energy excitations in finite and a pri- 
ori continuous systems of interacting electrons at low densities we have derived a lattice 
formulation for this problem. Non-orthogonal single electron states centred near the po- 
sitions where the electrons would be in the classical minimum energy configurations are 
orthogonalised, leading to contributions from double occupancies ~ UiU^riii in an effective 
Hubbard model description. Subsequent elimination of high energy states ( ~ Ui or ~ V%j , 
from nearest neighbor occupations) yields effective antiferromagnetic Heisenberg (eq. ( |3.30| ) 
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for v = 1 ) or charge-spin Hamiltonians (eq. (|3.28|) for v > 1 ) depending on the number v of 
classical minimum energy configurations. An extended version of the Hubbard model must 
be considered if highly symmetric geometries or ring processes are important. A crucial 
enhancement has been to enable "ring exchange" processes to be included, corresponding 
to the cyclic permutation of N > 2 electrons. Intermediate lattice points were introduced 
and the additional states corresponding to their single occupation were also eliminated by 
higher order perturbation theory, even though they are much lower in energy than states 
corresponding to double occupation. This yields generalized effective spin- 1/2 models con- 
taining up to (iV — l)-fold products of pairs of spin operators (see Eqn. |4.7| ). The magnitudes 
of the parameters for the superexchanges J and the ring exchanges K can be estimated 
semiclassically. 

Apart from providing a quantitative description of the low-energy spectra in quantum 
dots the mapping onto lattice models is, in itself, very illuminating. For example, in cases 
that can be mapped onto the antiferromagnetic Heisenberg model for the low energy states, 
the spectrum consists of a level multiplet of 'band width' NJ where J , the nearest 
neighbor superexchange integral, depends exponentially on the mean electron distance r s . 
The highest energy state within this multiplet is the spin polarized state, in agreement with 
what has been derived within the pocket-state approximation using less obvious arguments. 
All spectra obtained by the lattice descriptions coincide exactly with those found in the 



pocket-state |L6[ approximation and are in very convincing agreement with results gained 
by numerical diagonalisations [|U|,[JIJ. The computational effort required for the former is, 
however, considerably less compared to both of the other methods. This will enable much 
larger systems with iV > 20 to be investigated. 

In certain geometries of two dimensional quantum dots, v > 1 classical configurations 
of lowest electrostatic energy may exist leading to an effective spin-charge model instead of 
a spin only effective model. Whilst such situations are straightforward to deal with by the 
lattice method, we should point out that in real experimental situations the high symmetry 
will be lowered by a polarizable environment for the following reason. The energetically 
most favourable locations for the quantum dot electrons depend on the distribution of sur- 
rounding (non-conducting) charges, which themselves are influenced by the distribution of 
the (Wigner localised but under transport conditions conducting) dot charges. Selfconsis- 
tency allows the dot electrons to lower the total energy by adjusting their environment. For 
example, two electrons in a square will easily break the two- fold symmetry of the classical 
ground-states by polarizing their surroundings, leaving a diamond shaped configuration for 
the potential. The interplay between the surroundings and the granular electron density of 
the dot finally tends to lower the number of classical low energy configurations. Only config- 
urations with energies smaller than the band width ~ NJ eventually need to be considered 
explicitly. This plays a particular role e.g. in double-dot systems [[36|. Additionally, one 
should keep in mind that the shape of the quantum dot also depends on voltages applied to 
side gates |p7[| , so that the fine structure spectra may change with gate voltage. 

The lattice description simplifies considerably the determination of excitations in quan- 
tum dots. These excited levels may be deduced implicitly from non-linear transport ex- 
periments [H. Furthermore, transport is qualitatively influenced by the total spins of the 
many-electron states by two types of spin blockades P,p8|. 

The occurrence of a non-minimal total spin for the ground state in a finite electron 
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system is a direct consequence of the geometry of the boundary. For systems which can 
be represented by an antiferromagnetic Heisenberg model, this relationship is obvious (c.f. 
discussion of the square dot with N = 5 in section [V]). Apart from the electron number the 
spin values of ground and of excited states should depend very sensitively on the shape of 
the quantum dot. In future research this will be investigated systematically to allow specific 
suggestions for experimental design. 

Another highly relevant topic for real experimental situations is to apply the mapping 
derived in the present work to study the interplay between mutual interactions of the quan- 
tum dot electrons and impurities. No reliable picture exists presently, which demonstrates 
our lack of understanding for the persistent currents observed in mesoscopic rings. Theory 
is still not capable of explaining even the order of magnitude of the high values found exper- 
imentally. It is known that one dimensional models of rings are insufficient to describe this 
problem [p9j] . Furthermore, evidence exists that the electron spin is important 0,^]. Our 
mapping onto lattice models demonstrates that calculations based on the Hubbard model 
pdH43| are not significantly influenced by the lattice but indeed reproduce the low energy 
properties of the original continuous problems in presence of strong interactions correctly if 
the parameters, particularly the filling of the former, are interpreted accordingly. Advantage 
can be taken of the comparably large electron numbers that can be dealt with by our method, 
which would be otherwise intractable. This will allow the investigation of multiple-channels 
in two-dimensional rings of finite widths in the presence of impurities and flux. 
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TABLES 



TV 


5 


m -^Ground state 


2 
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J(2) 


3 


1/2 





3 


1/2 


J(3) 


3 


3/2 


(3/2)J< 3 > 


4 








4 


1 


(1 + y/3- ^2)J (4) /2 


4 


1 


(1 + v/3) JW/2 


4 





(2^3)J (4) /2 


4 


1 


(1 + V3 + >/2)J< 4 >/2 


4 


2 


(3 + ^3)J (4) /2 



TABLE I. Analytical values for the fine structure spectrum E^p of N interacting electrons 
in a one-dimensional hard wall box within PSA for N < 4 . S is the total spin of N Fermions 
with s = 1/2 . The excitation energies — E^l und statc , are given in units of J (iV) . 
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FIGURES 



FIG. 1. Classical low-energy configurations for two electrons in a 2D square quantum well. The 
first two configurations are lowest in energy whilst the remaining four are V higher in energy due 
to Coulomb repulsion. All other configurations have higher electrostatic energy. 

FIG. 2. (a) Classical ground-state configuration for three electrons in a triangular quantum dot. 
(b) Classical path for the exchange of two electrons, (c) Classical path for the cyclic permutation 
of all three electrons. Particles are shown in intermediate positions along their trajectories. 

FIG. 3. Lattice points for the triangular quantum dot. In the ground manifold all three elec- 
trons will be close to the positions indicated by the circles. Excited (intermediate) states correspond 
to one or more electrons in localised orbitals centered on the crosses. 

FIG. 4. Superexchange processes, (a) Fourth-order (~ t 4 /V 2 U). (b) Fifth-order (~ t 5 /V 4 ). 

FIG. 5. A six-step ring process for three electrons in a triangular quantum dot. The elec- 
trons would move continuously and simultaneously in the direction of the arrows in the WKB 
approximation . 

FIG. 6. Fine structure multiplets of a quasi one-dimensional quantum dot, for N = 3, . . . , 6 
as obtained directly from the pocket-state approximation and in exact agreement with the eigen- 
solutions of the effective spin Hamiltonian. The dashed lines were obtained by direct numerical 
diagonalisation of the interacting electron problem. The TV-dependence of j( N ) is not considered 
and J has been adjusted to normalize the 'bandwidth' of the multiplets. 

FIG. 7. Fine structure spectra of (a) N = 3 , (b) N = 4 and (c) N = 5 electrons in a 2D 
square as obtained within the pocket-state approximation. In (a) the dominant tunneling process 
(t) corresponds to nearest-neighbor hopping of an electron whereas in (b) and (c) the dominant 
tunneling process is nearest-neighbor exchange (J). All spectra are identical with those obtained 
from the corresponding charge-spin Hamiltonian. 
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